PASJ: Publ. Astron. Soc. Japan , 1-??, 
© 2008. Astronomical Society of Japan. 

Deep 15 fim AKARI observations in the CDFS : estimating dust 
luminosities for a MIR-selected sample and for Lyman Break Galaxies 
and the evolution of Ldust/ Luv with the redshift. 

Denis Burgarella * 

Laboratoire d'Astrophysique de Marseille, Observatoire Astronomique Marseille Provence, Universite dAix-Marseille, CNRS 

denis. burgarella@oamp.fr 
Vcronique Buat 

Laboratoire dAstrophysique de Marseille, Observatoire Astronomique Marseille Provence, Universite dAix-Marseille, CNRS 

veronique. buat@oamp.fr 
Tsutomu T. Takeuchi 

Institute for Advanced Research, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan 

takeuchi@iar. nagoya-u. ac.jp 
Takchiko Wada 
Department of Infrared Astronomy, ISAS, JAXA 
wada@ir. isas.jaxa.jp 
and 

Chris Pearson 

Space Science and Technology Department, CCLRC Rutherford Appleton Laboratory, Oxford 
Department of Physics, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta T1J 1B1, Canada 

c.pearsonQrl. ac. uk 

(Received 2008 May 23; accepted 2008 November 6) 
Abstract 

Deep observations of the Chandra Deep Field South have been secured at 15/im with AKARI/IRC 
infrared space telescope. From these observations, we define a sample of mid infrared-selected galaxies 
at 15fim and we also obtain 15/im flux densities for a sample of Lyman Break Galaxies at z ~ 1 already 
observed at 24/im with Spitzer/MIPS. Number counts for the mid infrared-selected sample show a bump 
around a 15/zm flux density of 0.2m Jy that can be attributed to galaxies at z > 0.4 and at z > 0.8 for 
the fainter part of the bump. This bump seems to be shifted as compared to other works and a possible 
origin can be the Cosmic variance. Thanks to this dataset, we have tested, on the two above samples 
z ~ 1, the validity of the conversions from monochromatic luminosities vj v at a rest- frame wavelength of 
8fim by a comparison with total dust luminosities estimated from Spitzer rest-frame 12/^m data that we 
use as a reference. We find that the 8fim dust luminosities are not all consistent and that some of them 
are better when compared to Ldust evaluated from longer wavelength luminosities. We also find that the 
rest- frame 8//m luminosities provide globally good estimates of L^ust- By comparing our data for the two 
samples to several libraries of spectral energy distributions, we find that models can explain the diversity 
of the observed /24//15 ratio quite reasonably for the MIR-sclcctcd sample and better for the LBG sample 
which are less dispersed than the MIR selection. However, when we analyse the luminosity dependence 
of this ratio, we find important discrepancies. Finally, we revisit the evolution of L^ustj Luv ratio with 
the redshift z by re-calibrating previous Ldust at z ~ 2 based on our results and added new data points at 
higher rcdshifts. The decreasing trend is amplified as compared to the previous estimate. 

Key words: galaxies: formation - galaxies: evolution - galaxies: high redshift - infrared: galaxies - 
cosmology: observations 



1. Introduction 

The star formation rate (SFR) of galaxies can be seen 
as some kind of Sesame that opens up the door of a quan- 
titative understanding of the formation and the evolution 
of galaxies at all rcdshifts. It is one of the main tools of 
modern cxtragalactic astrophysics and astronomers world- 
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wide have been looking for the best one(s) for several 
decades. 

All young stars in galaxies emit ultraviolet (UV) pho- 
tons that can be collected by telescopes to estimate the 
SFR of these galaxies (e.g. Lcithcrcr ct al. 1999). 
However young stars might be embedded in dust that 
prevents part of the UV photons to reach our telescopes. 
These UV photons are scattered and absorbed by dust 
grains and reprocessed into infrared (IR) photons (e.g. 
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Buat & Xu 1996) which poses a problem to obtain a 
reliable SFR than can be seen as a total one. 

Several methods have been proposed to evaluate the to- 
tal amount of SFR, including that part hidden in dust 
which cannot be measured by UV/optical telescopes. 
Meurer et al. (1999) proposed a version based on a cor- 
rection for the dust attenuation that would be applied to 
the UV flux. This correction relies only on UV data and 
is estimated from the UV slope f3 (assuming fxocX 13 ): the 
redder the UV slope the larger the correction. However, 
this method seems not to be fully reliable for any type of 
galaxies (Buat ct al. 2005) or, at least, provides an esti- 
mate that is not very accurate. Usually, the calibration is 
performed from infrared far-IR (FIR) data but this hap- 
pens not to be that simple because of various technological 
challenges: the low transparency of the atmosphere that 
implies to launch space telescopesl to observe above a few 
micrometers, the necessity to develop efficient detectors 
in this wavelength range and, of course, the mandatory 
cooling of these telescopes. Moving from the FIR to the 
mid-IR (MIR) looks like a good idea since the MIR emis- 
sion is easier to measure but, again, there might be issues 
that need to be carefully checked. First, the MIR range 
can be subject to a contamination from stellar emission 
(especially when observing high rcdshift objects). Next, 
some polycyclic aromatic hydrocarbons (PAHs) emission 
features are present in this range that are correlated to 
the metallicity (Engelbracht et al. 2008). Then, most 
of the SFR relationships have been calibrated from (or 
at least mostly rely on) observations of local galaxies. We 
can wonder whether they are also valid at higher redshifts. 
Indeed, the spectral energy distributions (SEDs) of galax- 
ies can evolve with the redshift and we do not definitely 
know, for instance, whether or not low and high redshifts 
share the same PAH emissions in the MIR. 

Moreover, even when higher redshift objects have been 
accounted for (e.g. Bavouzet et al. 2008), these galax- 
ies are not UV-selected galaxies. So, we have to check 
whether the relationships dcvcloppcd from IR samples are 
valid for UV-selected samples. 

Our objective, in this paper, is to use two samples of 
galaxies: one selected from deep 15/im observations of 
the Chandra Deep Field South (CDFS) and one sample 
of Lyman Break Galaxies (LBGs) at z ~ 1 from GALEX 
(Burgarella et al. 2006 and Burgarella et al. 2007) 
to compare estimates of the dust luminosities (computed 
from rest-frame 8/im - i.e. AKARI observations at 15/zm 
- using formulae provided by several authors) to the ones 
evaluated from Spitzer observed flux densities at 24/xm 
(i.e. rest- frame 12/im at z = 1). Then, we will test the 
MIR SEDs through a comparison of the observed /24//15 
flux ratio with several libraries of SEDs. 

Recent works by Reddy et al. (2006) and Burgarella 
et al. (2007) suggest that the dust-to-UV luminosity 
ratio (Ldust/Ljjv) would decrease from z = to z = 2 for 
a galaxy with a given luminosity L^ust- However, this 
effect is based on a comparison of datasets at different 
redshifts and we need to revisit this effect after checking 
the conversions from MIR luminosites to Ldust- 



In Sect. 2, we will present the observations of the 
Chandra Deep Field South with AKARI/IRC camera 
through the L15 filter centered at 15/im. Sect. 3 will 
present the two samples that we will use in this paper. 
Then, we will test several calibrations of 8/xm MIR lu- 
minosities to total dust luminosities. In Sect. 4, we will 
study the evolution of IR-to-UV luminosity ratio with the 
rcdshift by accounting for the results found in the paper. 
Finally, Sect. 5 will present the spectral energy distribu- 
tion libraries that we will use. We will try to understand 
better the physics of the two galaxy samples by analysing 
the /24//15 ratio by a comparison with the models and 
templates. 

We assume a cosmology with Hq = 70 km.s^ 1 .Mpc^ 1 , 
Q M = 0.3 and Ov = 0.7. 

2. Observations and Data Reduction 

The AKARI satellite (Murakami et al. 2007) is a 
Japanese space mission dedicated to infrared astrophysics 
and was launched on board JAXA's M-V8 Launch Vehicle 
on February 22, 2006 Japan Standard Time, JST. AKARI 
has a 68.5 cm cooled telescope with two focal plane instru- 
ments, the Far-Infrared Surveyor (FIS) (Kawada et al. 
2007) and the Infrared Camera (IRC). The IRC consists 
of three cameras, the IRC-NIR, MIR-S & MIR-L covering 
1.7-26 fim in 9 bands with fields of view (FoV) of 10' x 10' 
(Onaka et al. 2007). Although the NIR/MIR-S cameras 
share the same field of view, the MIR-L camera is sepa- 
rated on the sky by 20 arcmin. AKARI observations are 
segregated into Large Surveys (LS) and guaranteed time 
Mission Programs. In addition, 10% of the total observa- 
tion time was reserved for Open Time observations from 
ESA member institutes. In this work, we report on our 
Open Time program (IRLBG, P.I. D. Burgarella) to im- 
age the Chandra Deep Field South (CDFS: RA=53.0465, 
Dcc=-27.8307, J2000) in the AKARI-IRC 15/im band 
(the L15 band covering the wavelength range 12.6 - 
19.4/im). A total of 5 pointed observations were taken 
on the 31st January and 1st February 2007, using the 
IRC Astronomical Observing Template (AOT) designed 
for deep observations (IRC05) in the configuration b;L, 
corresponding to the L15 filter on the IRC-MIR-L camera 
(See Table 1). In this configuration, each pointing consists 
of approximately 10 exposure cycles. Each exposure cycle 
consists of a single short (0.6s) exposure frame and three 
long (19s) exposure frames. For deep imaging, the short 
exposure frames are insignificant and are often discarded. 
Note that, for a target centered on the FoV of the MIR-L, 
data is also taken in the NIR/MIR-S cameras, albeit offset 
by ~20 arcmin with the FoV rotating around the MIR-L 
target position. The IRC05 AOT has no explicit dither- 
ing built into the AOT, therefore dithering is achieved by 
carrying out at least three pointed observations with this 
AOT at any given position. The dithering is achieved from 
the natural slight offset movement and slight rotation of 
the field of view between any pair of pointed observations 

The data were reduced by the standard IRC imaging 
pipeline (version 070911; see IRC Data User Manual by 
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Lorente et al. 2007). The pipeline corrects basic instru- 
mental effects, performing dark subtraction, linearity cor- 
rection, distortion correction and flat fielding. The super- 
dark frame version 061109 and super flat version 070912 
were used. In order to reduce the effect from the numerous 
hot pixels in the MIR-L detectors, hot pixels were identi- 
fied and set to a very large numeric value such that they 
would be identified and removed as cosmic ray events in 
a later processing stage of the pipeline. The pixel values 
were then interpolated using the surrounding pixel values 
after the flat fielding. 

AKARI observations suffer from the effects of stray 
light from Earthshinc and we have rejected any frames 
whose sky value had a measured standard deviation sig- 
nificantly larger than the mean between frames. Before 
combining the image from individual exposure frames for 
a single pointing, the frames have a median sky sub- 
tracted. Alignment of the frames for co-addition is usu- 
ally determined by the position of bright point sources 
in each frame. However, in case of the MIR-L chan- 
nel, the number of the point sources which had suffi- 
cient signal-to-noise were too few. Therefore, in order to 
align the MIR-L frames, we identified common objects in 
our MIR-L frame and in the Spitzcr 24/im FIDEL image 
(http://ssc.spitzer.caltech.edu/legacy/). 

Spurious events, such as cosmic rays and electronic 
ghost images are removed at the co-addition stage using a 
3 sigma clipping technique. For our data, given the large 
number of frames available for co-addition, we used the 
cedclip option in the pipeline, where the detector noise 
performance is assumed as opposed to the default option 
sigclip where the noise performance is estimated from the 
data itself. The average rather than the median value 
among the frames, for each pixel was used in order to 
further improve the signal-to-noise of the final image. For 
every pointing a total of 30 exposure frames were co-added 
corresponding to a total exposure time of 570s. The five 
individual pointings were then combined to make the final 
image (Fig. 1). 

3. Definition of the MIR-selected and the Lyman 
Break Galaxy samples 

The above AKARI image has been used to build two 
galaxy samples : a 15fj,m - selected (or cquivalently MIR- 
selected) sample defined from the photometry carried out 
on the image. We also measured the 15/im flux density 
of a sample of Lyman Break Galaxy (LBG) already de- 
fined from GALEX far ultraviolet (FUV) and near ultra- 
violet (NUV) images (Burgarella ct al. 2007). 

3.1. The MIR - selected sample 

The IRAF task DAOFIND has been used to identify 
a set of objects at a 3 — a level above the background. 
To minimise confusion errors, we have integrated over an 
aperture of 2 pixels using the IRAF task PHOT and then, 
this value has been extrapolated (by a factor of 1.645 es- 
timated from a PSF built from several point sources in 
the field) to an integrated flux density over 7.5 pixels to 



match the recommendation of the AKARI / IRC man- 
ual by Lorente et al. (2007). Then, the calibration pro- 
vided by (Lorente et al. 2007, i.e. /„ = fr.5 P ixels x 1-48) 
has been applied to the data to get flux densities in fiJy. 
Uncertainties on the flux densities are provided by PHOT. 
The estimate of the uncertainties also account for the 2.82 
% error on the calibration constant quoted by Lorente et 
al. (2007). A list of 323 MIR - selected objects is created. 
Then, we have cross-correlated this sample with COMBO 
17 (Wolf et al. 2004). Whenever several counterparts are 
found within the 2-arcsec radius, we select the best one 
by eye. The final sample that we use in this paper con- 
tains 239 objects that we also have cross-correlated with 
the 24/xm-selected sample from Buat et al. 2008. The 
dataset is presented Table 2. 

The redshift distribution of this sample is shown in Fig. 
2. It shows that galaxies up to z = 1.3 have been detected 
in the AKARI image. We select only objects at z < 1.2 
for this analysis because the quality of the highest redshift 
objects is too low. 

3. 2. The z ~ 1 Lyman Break Galaxy sample 

Burgarella et al. (2007) defined a sample of 420 LBGs 
at0.9<z<1.3 from ultraviolet images in the deep FUV 
and NUV exposures in the CDFS. In this paper, the au- 
thors find that this LBG sample shares its main char- 
acteristics with UV-selected samples. The LBG sample is 
formed of 80% not detected at 24/j,m (we called them Blue 
LBG = BLBGs) and 20% detected at 24^ra (Red LBGs = 
RLBGs). Note that this classification is only representa- 
tive of a difference in dust attenuation but the two of them 
share the same other properties in average (star formation 
history, morphology, etc.). However, those 20% of RLBGs 
contain 80% of the SFR density of LBGs. Burgarella ct 
al. (2007) also claim that the UV - optical spectral energy 
distributions (SEDs) of the z ~ 1 sample is very simular 
to the z ~ 3 LBG sample selected by Forster-Schreiber et 
al. (2004). 

One of the initial motivation of the AKARI program 
is to check whether we could use the rest-frame ~ 8{im 
flux density to estimate the hidden SFR of LBGs (or UV- 
selected objects) at high redshift. Several authors have 
calibrated this rest-frame ~ 8/xm flux density into a hidden 
star formation rate by using z = galaxies (e.g. Chary & 
Elbaz 2001, Takeuchi et al. 2005b, Reddy et al. 2006). 

We have used the coordinates from the z ~ 1 LBG sam- 
ple (Burgarella et al. 2007) as inputs for the IRAF 
task phot, in the same way as for the MIR-selected sam- 
ple. After cleaning up the results from the photometry, 
27 LBGs are found to have a Ibfim flux density within 
the range 17.39 < s L5Aim \[iJy\ < 227.78. Out of these 28 
LBGs, 9 have been detected at 24{im in the redshift range 
0.9 < z< 1.2 used in the rest of the paper. This small LBG 
dataset only contain RLBGs which are IR-bright objects. 
Moreover, since these galaxies are LBGs, they have been 
detected in the GALEX near-UV bands. This insures that 
we are not using higher redshift galaxies as z ~ 1 ones : at 
z ~ 1.4, the Lyman break would be beyond the red limit 
of the near-UV filter and no flux at all should be detected 
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Table 1. Observation Log 



Target id AOT OBS-DATE RA DEC 

4040001-001 IRC05 b;L 2007-01-31 04:05:30 53.0465 -27.8307 

4040001-002 IRC05 b;L 2007-01-31 05:44:49 53.0465 -27.8307 

4040001-003 IRC05 b;L 2007-01-31 23:57:15 53.0465 -27.8307 

4040001-004 IRC05 b;L 2007-02-01 01:36:34 53.0465 -27.8307 

4040001-005 IRC05 b:L 2007-02-01 03:15:52 53.0465 -27.8307 
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Fig. 1. Combined AKARI / IRC 15 fim frame 



Table 2. 15/^m-selected sample 



COMBO 17 Id RA Dec z Si 5Mm errs 15flm s 2 4 M m errs 2 A^m L dust 



Table 3. z ~ 1 Lyman Break Galaxy sample 



COMBO 17 Id RA Dec z sis^m errsi 5flm s 2 4 M m errs24 M m £d«st 
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Fig. 2. The redshift distribution is evaluated by using COMBO 17 redshifts. 



in this filter. The dataset is presented in Table 3. 
4. The 15 fim number counts 

In this section, we show the number counts of AKARI 
15-fim sample in the CDFS. The total number of the 15- 
fim galaxies in this field is 239. Cumulative number counts 
(NC) of these galaxies are shown in Fig. 3 (dotted irregular 
step function). The source density is converted to the 
units of [sr _1 ] in Fig. 3. 

We can see the effect of a possible incompleteness below 
a flux density of ~ 100 fiJy. We have examined the incom- 
pleteness of the source extraction by adding additional 
sources with various flux densities on the original observed 
AKARI 15-/mi image (IRAF task ADDSTAR) and run- 
ning the source extracting task (IRAF DAOFIND). We 
can estimate the incompleteness by measuring how much 
fraction of the added artificial sources are detected by the 
source extraction procedure. We show the result of this 
examination in Fig. 4. The open squares with error bars 
are the measures rate of retrieving artificial sources. To 
have a continuous function of the completeness or "miss- 
ing rate" , we fit the measured value with a function of the 
Fermi-Dirac form, 



1 



exp 



(1) 



where we found that Sc = 57.3 and <fi = 21.1. By the 
intersection of this function and 0.8 (horizontal line in 
Fig. 4), we found that the 80-% flux completeness limit 
is 85 /iJy. With eq. 1, we can correct the raw NC and 
obtain the completeness-corrected NC. The blue regular 
step function in Fig. 3 is the corrected NC. 

Based on the integrated NC, we obtain the differen- 
tial NC in a standard manner. We show the Euclidean- 
normalized differential NC in Fig. 5, i.e., dN/dSi§ mul- 
tiplied by S15 2 ' 5 . The open triangles represent the raw 
differential NC from our CDFS dataset and red diamonds 
represent results from Wada et al. (2007) in the NEP. 
The errors are Poissonian. Again, with eq. 1, we can cor- 
rect the incompleteness effect. The open squares in Fig. 5 
show the corrected differential NC. Error bars are cal- 
culated assuming Poissonian fluctuations in the original 
counts, and multiplied the incompleteness for those with 
corrected counts. The vertical dotted line is the same as 
Fig. 3. As we can see, both differential NC are consistent 
but show a variation around the peak of CDFS observa- 
tions and the AKARI Performance Verification (PV) one. 
However, the cosmic variance in the CDFS (GOODS, to 
be precise) is estimated to be at least 60 % in density- 
related measurements (Somerville et al. 2004). Since the 
area is smaller than the whole GOODS area, we should 
have larger variance both in CDFS 15^m and in the PV 
field. In addtiion, is known that the CDFS hosts clusters 
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Fig. 3. The integral number counts (NC) of galaxies in the CDFS observed with AKARI IRC MIR L15 band. The green irregular 
step function (dotted line) shows the raw NC of the galaxies, while the regular step function (solid line) represents the NC corrected 
for the completeness. The vertical dotted line show the 80-% flux completeness limit, which is found to be 85 /^Jy (see text). 
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Fig. 4. Completeness of the source extraction of our CDFS AKARI 15-/im image. The open squares with error bars are the measures 
rate of retrieving artificial sources, and the solid curve shows the fitted function. We found that the 80-% flux completeness limit is 
85 /Jy. 
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at z ~ 0.7 (e.g. Adami et al. (2005)). So, the observed 
discrepancy would not be statistically significant with re- 
spect to the cosmic variance (especially at 15/im i.e. 7{im 
rest-frame at z <~ 1 which is sensitive to any star formation 
in the clusters) and to other datasets presented in Figures 
12 and 13 from Wada et al. (2007). Quite surely we are 
observing the same population of galaxies as they do. 

The most prominent feature of the differential NC is the 
bump seen around a flux density of 200 yuJy. In order to 
see what population consists of this bump, we divide the 
differential counts into galaxies with three redshfit bins: 
< z < 0.4, 0.4 < z < 0.8, and 0.8 < z. We show the sub- 
divided differential counts in Fig. 6. The smallest open 
squares, intermediate ones, and the largest ones represent 
galaxies at < z < 0.4, 0.4 < z < 0.8, and z > 0.8, respec- 
tively. The broken solid line is the total differential NC, 
the same as shown in Fig. 5. It is clear that the bump 
is dominated by galaxies with z > 0.4. Especially, the 
fainter part of the bump consists of galaxies at z > 0.8. 
This may mean that galaxy population observed at 15 /im 
experienced a significant evolution, as suggested by previ- 
ous studies based on number counts (Pearson & Rowan- 
Robinson 1996; Pearson 2005; Takeuchi et al. 2001; Xu 
et al. 2001). 

5. SED Libraries 

5.1. Spectra, templates and models in the Mid Infrared 

We have used several libraries of templates to be com- 
pared to the data. First, three SED libraries commonly 
used to analyse MIR properties of local and distant galax- 
ies: those of Chary & Elbaz (2001) (CE), Dale & Helou 
(2002) (DH) and Lagache et al. (2004) (LDP). These dif- 
ferents sets have been developped and optimized for differ- 
ent purposes: the CE library reproduces the local correla- 
tions between different luminosities from 7 to 850 /im, the 
DH templates describe the infrared SEDs of local normal 
star forming galaxies as a function of a single parameter, 
fool f wo, and the LDP templates of starburst and normal 
galaxies are built and optimized to reproduce statistical 
quantities like number counts, rcdshift distributions and 
the cosmic infrared background. The original templates 
for starburst galaxies (Lagache et al. 2003) have been 
modified in the MIR to reproduce the new Spitzer num- 
ber counts with a decrease of the amplitude of the 11 
/zm feature by a factor ~ 2 (Lagache et al. 2004). The 
LDP templates for normal galaxies all exhibit the same 
shape and are only scaled in luminosity. They are similar 
to DH templates with intermediate feoff wo values and 
thus, moderate luminosities. So, we do not reproduce all 
of them here. 

We have also considered recent templates built from 
the low-resolution Spitzer IRS spectra. From their anal- 
ysis of IRS spectra for 59 nearby galaxies Smith et al. 
(2007) have proposed five templates from 5 to 38 fim, one 
is a luminosity weighted average of all spectra, the other 
ones represent averages of spectra with different strengths 
of the features inside the PAH bands (weak AGNs are 
also included in the sample). To enlarge the range of 



observed spectra we have also included in our set of tem- 
plates the nine classes of spectra proposed by Spoon et 
al. (2007). Their analysis is based on the observation 
of a sample of AGNs, starburst galaxies together with 
LIRGs and ULIRGs. Their MIR templates exhibit very 
different spectral shapes from continuum dominated spec- 
tra for AGNs to spectra dominated by the 10 fim silicate 
trough and PAH-dominatcd starburst spectra. Last but 
not least, we have also included the starburst template 
of Brandl et al. (2006) obtained by averaging the IRS 
spectra of 13 nearby starburst galaxies. 

The different sets of templates are illustrated in Fig. 7 
and Fig. 8 where a few representative spectra for each li- 
brary are plotted. For the CE and LDP libraries (Fig. 7) 
we have chosen the templates corresponding to galaxies 
with L dust equal to 10 11 and 10 12 L©. The DH tem- 
plates are given as a function of feoff wo- F° r the pur- 
pose of our comparison between the differents libraries, 
it is convenient to attribute to them a total dust lumi- 
nosity. At this aim, we follow Marcillac et al. (2006) 
who fitted a linear relation between the logarithms of 
Ldust and feoff wo valid for the nearby universe. We 
have plotted the templates with log(/6o//ioo)=-0.2 and 
l°g(/6o//ioo)=-0.1 which correspond respectively to a to- 
tal dust luminosity of the order of 10 11 and 10 12 LQ. The 
templates obtained from Spitzcr/IRS spectra are plotted 
in Fig. 8: the luminosity- weighted average spectrum and 
one other spectrum (template 3) from Smith et al. (2007), 
the templates 1A (continuum dominated), 3B (prominent 
silicate absorption), 1C and 2C (PAH dominated) from 
Spoon et al. (2007) and the starburst template of Brandl 
et al. (2006). 

5. 2. The rest-frame 8fim emission redshifted at z = 1 and 
z = 2 

Several paper provide a calibration of the rest-frame 
8/im luminosity into a dust luminosity (e.g. Chary & 
Elbaz 2001, Takeuchi et al. 2005a, Reddy et al. 
2006). Given the availability of deep 24/ito observations, 
these data have been used to estimate dust luminosities 
at high redshifts (e.g. Reddy et al. 2006). However, 
if we wish to apply the z = conversions to higher red- 
shift galaxies, we must account for corrections because 
the Spitzer /IRAC 8/im filter at z = does not sample ex- 
actly the same part of the SED than the AKARI lhfim 
filter or the Spitzer/IRS \6fim filter at z = 1 or than the 
Spitzer/MIPS 24/im filter at z = 2 (Fig. 9). Moreover, 
these corrections will depend on the SED of the observed 
galaxies. Corrections are evaluated by convolving several 
templates and spectra from Chary & Elbaz (2001), Dale 
& Helou (2002), Lagache et al. (2004), Smith et al. 
(2007), Spoon et al. (2007) and Brandl et al. (2006). 
By using several templates and observed SEDs, we aim 
at taking into account the variety of observed SEDs. The 
observed Ldust distribution extends from 10 9 to 1O 115 L0, 
so only models in this luminosity range are used to build 
Fig. 10. The corresponding median corrections are listed 
in Tab. 4. We find that the dispersion in the correction 
factors ((p futer ) is relatively low (a few percents) for each 
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Fig. 5. Differential NC of the CDFS AKARI 15-/^m galaxies. The green open triangles represent the raw differential NC, while 
the blue open squares show the differential NC corrected for the incompleteness. Error bars are calculated assuming Poissonian 
fluctuation. The red open diamonds are from Wada ct al. (2007) in the NEP. 



filter and close to 4>futer = 1, which brings some confidence 
in using L^ust (8^m) . The correction factors will be used 
to compare the different estimates of Ldust at various red- 
shifts. The redshift range of our objects is 0.9 < z < 1.3. 
In this range, PAH features can move in and out of the 
filter. We need to check whether this can influence the 
results presented in this paper and, if so, what is the size 
of any effects. To perform this test, we have convolved 
SEDs from Dale & Helou (2002) and from Chary & Elbaz 
(2001) with the AKARI 15/oti filter. Since all filters will 
sample the same rest-frame wavelength range, this should 
address our concern. The values (s z= o.9 — s z =i.i)/s z =i.i 
and (s z= i.3 — s z= i.i)/s z= i.i are computed and shown Tab. 

5. SEDs from Chary & Elbaz (2001) provide somewhat 
larger relative change by about 16% as compared to the 
5% from Dale & Helou (2002) SEDs. We can conclude 
that the redshift range is small enough to keep within the 
uncertainties even by assuming a single redshift. 

6. Testing the rest-frame 8fim luminosity conver- 
sions into total dust luminosity 

Estimating the star formation rate of high redshift 
galaxies always is an issue, especially when one wishes 
to estimate the hidden star formation rate or, in other 
words the dust luminosity. Of course, if FIR data are 
available as for galaxies detected in IR surveys, the task 
is more straightforward since we are close to the maximum 
of the dust emission, but these data are, most often not 



observed for ultraviolet-selected galaxies like LBGs. Even 
though 8^m luminosities are calibrated into star forma- 
tion rates, we can wonder if we can apply these calibra- 
tions and, if so, what is the associated uncertainty. We 
directly compare in this paper calibrations for high red- 
shift galaxies which are likely to share common physical 
characteristics in terms of dust temperature distribution, 
metallicity, PAH contribution, etc. 

Caputi et al. 2007 performed such a comparison at 
z = 2 on relationships derived by convolving templates 
with the transmission of the Spitzer MIPS-24^m band. 
Daddi et al. (2007), Papovich et al. (2007) also per- 
formed a similar job at z = 2. Here, wc use a different 
approach : we use our LBG selected at 0.9 < z < 1.2 
for which we have rest-frame 8/im data (from AKARI 
IRC 15/iTO-band observations) and rest-frame I2(im data 
(from Spitzer MIPS-24^m band observations). These 
LBGs have 11 < LogLdust < 12 which means that they are 
Luminous IR galaxies (LIRGs with 11 < LogLdust < 12). 
A direct comparison of the two estimates allows us to 
quantify the quality and the uncertainty of the former, af- 
ter making the hypothesis that the latter provides a good 
estimate at z ~ 1. The conversion is largely dependent 
on the assumptions made or the templates used (e.g. Le 
Floc'h et al. 2005). By using directly observed data, we 
wish to verify what the validity of Ldusti^fxm) is to esti- 
mate Ldust on actual high redshift galaxies at z = 1 . This 
work therefore bears some similarities with Marcillac et 
al. 2006 at z = 1 but our samples are difference and we 
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Fig. 6. Differential NC of the CDFS AKARI 15-/xm galaxies, subdivided into galaxies with three redshift bins (completeness 
corrected). The smallest open squares represent the galaxies at < z < 0.4, the intermediate ones are at 0.4 < z < 0.8, and the largest 
ones are 0.8 < z. The broken solid line is the total differential NC, the same as shown in Fig. 5. Clearly, the bump is dominated by 
galaxies with 0.4 < z. 



Table 4. To use rest-frame Sfim conversions from monochromatic luminosities to Ldusti we must apply some corrections <j> filter to 
the photometry. The corrections (f^, orr = f£ bs /<j> filter) f° r 2=1 for galaxies observed in the AKARI/IRC L15 15/im filter and in 
the Spitzer/IRS I6fim filter and for z = 2 for galaxies observed in the Spitzer/MIPS 24fim filter are listed here. 



Table 5. We have computed the ratio of the flux densities in the AKARI/IRC L15 15/xm filter using Dale & Helou (2002) and 
from Chary & Elbaz (2001) (s z =0.9 — Sz=i.i)/sz=i.i and (sj-i.3 — s z — i.i)/s z =i.i to test if the influence of PAH moving with the 
redshift of the observed galaxies might change the observed flux densities and produce an effect only due to K-corrections. 



Filter z = 1 AKARI 15/ito z — 1 Spitzer 16/ito z = 2 Spitzer 24/ito 
(j> filter 0.94 ±0.02 1.00 ±0.03 1.12 ±0.03 



SEDs 
Dale & Helou (2002) 
Chary fc Elbaz (2001) 



(s z =0.9 - S z =l.l)/S z= i.i (s z = 1.3 - Sz=l.l)/s z= l.l 



-0.092 ±0.047 -0.033 ±0.013 

-0.161 ±0.028 +0.016 ±0.035 
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Fig. 7. Sets of templates at L^ust = 10 11 LO (solid lines) and -L^ust = 10 12 L0 (dashed lines)) for the libraries of Chary & Elbaz 
(2001) (red), Dale & Helou (2002) (green) and Lagachc ct al. (2004) (blue), all the spectra are normalized at 10 /im, see text for 
details 



add the LBG specificity. 

We define Ldust as the total infrared luminosity from 
8/im to 1000/im. Another definition uses wavelengths in 
the range to 1000/ito. Bavouzet ct al. (2008) give 
the ratio £5-1000^/^8-1000^™ = 1.07±0.04. In this pa- 
per, we prefer using an empirical approach in this pa- 
per by using conversions that are mainly based observed 
datasets and do not rely on (or at least minimize the use 
of) templates or models. Theefore, to estimate Ldust from 
monochromatic MIR lumininositics, we compute two es- 
timates of Ldust at a rest-frame wavelength of Yl[im : 
Chary & Elbaz (2001) and Takcuchi et al. (2005a). By 
using these two formulae, we show that our results are not 
very dependent on the assumed 12^m conversion and we 
assume that Ldust (12/zm) can be used as a safe reference. 

Chary & Elbaz (2001) provide a formula estimated 
from a sample of local 300 IRAS Bright Galaxy Sample 
(Soifer et al. 1998). Almost all these galaxies are detected 
at all IRAS wavelengths from 12 to 100 [im, which means 
that they form a far-infrared selected sample dominated 
by relatively cool (i.e. Tdust ~ 50K) dust galaxies. 

Chary & Elbaz (2001) relationship corresponds to : 



T — fl QQ+0.38 7- 1.094 fn\ 
Ldust — U.89_ .27-^12^771 \ l ) 

Takcuchi et al. (2005a) also provide conversions from 
several mid- to far-infrared luminosities into total dust lu- 
minosities Ldust- The initial sample built by Takeuchi et 
al. (2005a) is extracted from the IRAS PSCz (Saunders 
et al. 2000) that contains 15411 galaxies with redshifts. 
Then, 1420 galaxies with good quality flux densities at all 
IRAS wavelengths were down-selected. Star forming qui- 
escent, normal galaxies dominate this sample. However, 
the authors found that their relationship provide a good 
estimate for a wide range of spectral energy distributions. 
The conversion is estimated by comparing monochromatic 
luminosities with Ldust computed using Sanders et al. 
(1996) calibration (the four IRAS wavelengths are used) 
where Ldust is computed by fitting a single temperature 
dust cmissivity model (with an emissivity e = to the 
flux in all four IRAS bands. 

Takcuchi et al. 2005a conversion (with a dispersion 
a = 0.154) follows the relationship : 

LogLdust = 1.02 + 0.972LogL 12fim (3) 

Then, for the conversion from 8fim monochromatic lu- 
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minosities to Ldust we made use of Reddy et al. (2006), 
Bavouzet et al. (2008) and Wu et al. (2005). 

Reddy et al. (2006) defined an optically-selected sam- 
ple with spectroscopic redshifts at z ~ 2. For these galax- 
ies, X-ray data are available through stacking and we have 
indications that the X-ray luminosities are correlated to 
the total far-infrared ones (e.g. Rosa-Gonzalez et al. 
2007) for galaxies with 8 < LogLdust [Lq] < 12. X-ray lu- 
minosities are proxies for the total infrared luminosities 
and thy provide a good check for the conversion that is 
not dependent on templates. To quantify the relationship 
between 8/zm luminosities and total infrared luminosities, 
Reddy et al. (2006) made use of Elbaz et al. (2002) data 
for a sample of 149 galaxies that contains ISO data and 
IRAS 60 and 100/im measurements. 

Reddy et al. 2006relationship is : 

W = (1.24±0.35)i 8Mm (4) 

Bavouzet et al. (2008) samples contain galaxies with 
spectroscopic redshifts and COMBO 17 photometric red- 
shifts observed at 8/ito, 24/ito, 70^m and 160/im by 
Spitzcr. This sample can be considered as a far-infrared 
selected sample with redshifts are up to z = 1 , but most of 
them have lower redshifts in the range z ~ 0.2 — 0.3. They 
have added 13 stacked data points that are similar to a 
mid-infrared selection and 3 stacked point from Zheng et 
al. (2007) that are representative of an optically-selected 



sample. Bavouzet et al. (2008) integrated the spectral 
energy distributions in rectangles to avoid biases to esti- 
mate total dust luminosities. 

Finally, the relationship finally defined by Bavouzet et 
al. (2008) (with an uncertainty of 37 %) is : 

Ldust = 450.9 x (ig^rn) ' 83 - (5) 

Wu et al. (2005) exploited the Spitzer First Look 
Survey imaged by Spitzcr IRAC and MIPS and also cov- 
ered by the Sloan Digital Sky survey (Stoughton et al. 
2002), which provides Ha fluxes for this galaxy sample. 
One of the advantges of having some spectroscopy is that 
narrow-lines AGNs can be removed from the sample. All 
the objects in this sample are star-forming galaxies. The 
Ha fluxes are corrected for dust attenuation computed 
from the Balmer decrement. A good correlation is ob- 
served between the 8/im luminosities and the Ha lumi- 
nosities. To deduce the relationship between v.fsum and 
Ldust, we convert their star formation rate into LogLdust 
with Kcnnicutt (1998) SFR formula. 

Wu et al. (2005) relationship corresponds to : 

Ldust = 0.566L 8Mm (6) 

The two left panels of Figs. 11 and 12 refer to the LBG 
sample while the two right panels are related to the CDS 
AKARI MIR-selected sample. 

The first obvious point is that all conversions presented 
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in Figs. 11 and 12 are not fully consistent and we can see 
differences on the resulting 8^m conversions into Ldust- 
We clearly observe a stratification of the 8^tm conversions 
(as compared to the 12/im ones : Chary & Elbaz 2001 and 
Takeuchi et al. 2005a) from Reddy et al. (2006) that 
over estimate Ldust both for LBGs and for the AKARI 
sample to Bavouzet et al. (2008) conversion that pro- 
vides globally the best estimate and to Wu et al. (2005) 
conversion that slighly under estimates Ldust- The con- 
version formula from Wu et al. (2005) is calibrated from 
Ha data corrected for dust attenuation using the Balmer 
decrement. Obtaining an under estimated Ldust is not un- 
expected since Ha observations might miss some deeply 
enshrouded emissions. 

We see a difference between the LBG and the MIR- 
selected diagrams. While for the former, the slopes of all 
the conversions are close to 1 (see Tab. 6), they are shal- 
lower for the MIR selection. It is interesting to notice that 
the relationship is almost perfect between Chary & Elbaz 
(2001) and Bavouzet et al. (2008). Both of them provide 
the same Ldust within the uncertainties. On the other 
hand, the slope is shallower for the MIR-selected sample. 
The origin of this difference seems to be due to the ex- 
tremes at low {Ldust ~ 10 11 !/©) and high (Ldust > 10 12 £q, 
i.e. ULIRGs) dust luminosities. Taking off high-Ldust 
galaxies in the fit does not change the slope which remains 



significantly different from 1. To reach similar slopes, we 
have to also take off the lowest luminosity galaxies from 
the MIR sample. Part of the difference in the slope can 
be explained by the fact that LBGs at Ldust ~ 10 n -L© 
have lower Ldust estimated from 8//m than the AKARI 
sample. A detailed analysis and a comparison to models 
will be performed in the next section devoted to the MIR 
SEDs and to the interpretation of the /24//15 flux ratio 
that might provide some explanation for this difference. 

7. The evolution of the Ldust/ Ljjv ratio of 
ultraviolet-selected galaxies with the redshift 

Burgarella et al. (2007) and Reddy et al. (2006) found 
a decrease of the Ldust/ Luv ratio of ultraviolet-selected 
galaxies with the redshift from z = to z = 2. However, 
since we found in the previous section that Reddy et al. 
(2006) over estimated Ldust, we must revisit the evolution 
of Ldust/ Ljjv assuming another conversion. Note that 
with the updated calibration, the median and mean Ldust 
for Reddy et al.'s sample are respectively 1.4 lO 11 ^© and 
1.4 10 11 ±4.3 10 11 L©. We also find in the previous sec- 
tion the conversion given by Wu et al. (2005) provides 
under estimated Ldust-, probably because of some missing 
emission too deeply enshrouded in dust. We will therefore 
use Bavouzet et al. (2008) conversion instead of Reddy 
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Table 6. Statistical information on the regressions between the rest-frame and I2fim L^ uat estimates for a linear fitting 

logL dust (l2^m) = a + logL dugt (12)im). 



Selection 


Ldust{8[im) vs. 






a 




b 


Correlation coefficient 


LBG 


Chary & Elbaz 


(2001) vs. 


Wu et al. (2005) 


1 


183 ±0 


095 


-2.400 ±1.085 


0.957 


LBG 


Chary & Elbaz 


(2001) vs. 


Reddy et al. (2006) 


1 


183 ±0 


095 


-1.731 ±1.085 


0.957 


LBG 


Chary & Elbaz 


(2001) vs. 


Bavouzet et al. (2008) 





982 ±0 


079 


0.192 ±0.900 


0.957 


LBG 


Takeuchi et al. 


(2005a) vs 


Wu et al. (2005) 


1 


331 ±0 


107 


-3.818±1.199 


0.957 


LBG 


Takeuchi et al. 


(2005a) vs 


Reddy et al. (2006) 


1 


331 ±0 


107 


-3.149±1.199 


0.957 


LBG 


Takeuchi et al. 


(2005a) vs 


Bavouzet et al. (2008) 


1 


105 ±0 


089 


-0.985 ±0.995 


0.957 


15/iTO 


Chary & Elbaz 


(2001) vs. 


Wu et al. (2005) 





685 ±0 


093 


3.267 ±1.057 


0.696 


15/iTO 


Chary & Elbaz 


(2001) vs. 


Reddy et al. (2006) 





685 ±0 


093 


3.936 ±1.057 


0.696 


15/iTO 


Chary & Elbaz 


(2001) vs. 


Bavouzet et al. (2008) 





569 ±0 


077 


4.896 ±0.877 


0.696 


15/iTO 


Takeuchi et al. 


(2005a) vs 


Wu et al. (2005) 





771 ±0 


104 


2.445 ±1.167 


0.696 


15/iTO 


Takeuchi et al. 


(2005a) vs 


Reddy et al. (2006) 





771 ±0 


104 


3.115± 1.167 


0.696 


15/iTO 


Takeuchi et al. 


(2005a) vs 


Bavouzet et al. (2008) 





640 ±0 


086 


4.214±0.969 


0.696 
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Fig. 11. The two panels compare the relative performances of several conversions of L&^m into Ldust as compared to Ldust (12/im) 
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the comparison of Wu et al. (2005) 8fim conversion with Chary & Elbaz (2001). The left-hand panel (a) contain the LBG sample 
while the right-hand panel (b) contain the CDFS AKARI Sfim selection. 
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Fig. 12. As for Fig. 11 but for Takeuchi ct al. 2005a. 
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et al. (2006) in the following to compute a new esti- 
mate of Ldust/ Ljjv at z = 2. This new conversion means 
that Log Ldust decreases by —0.37 ±0.05 for ultraviolet- 
selected galaxies in the luminosity range corresponding to 
LIRGs. The z = l Ldust I ^uv was computed by Burgarella 
et al. (2007) by using Chary & Elbaz (2001) conversion 
and we are, therefore, in the optimum case where both 
estimates from and 12/im Ldust are fully consistent. 
Note that replacing Chary & Elbaz (2001) conversion by 
Takeuchi et al. (2005a) one would change Ldust estimated 
from the rest-frame 12/xm luminosityes by —0.2dex. Fig. 
13 updates the one presented in Burgarella et al. (2007). 
The change applied to Rcddy et al. (2006) at z = 2 
amplifies the decrease with the redshift. We have added 
one point at z ~ 3 from Adelberger et al. (2000) esti- 
mated from their Fig. 11 at Ldust/ Ljjv = 0.7±0.7. The 
dust luminosity of their sample is about 10 10 — 2.10 12 L Q . 
Although the Ldust I Ljjv was estimated from the UV slope 
P (Meurer et al. 1999), Burgarella et al. (2007) suggests 
that it provides a rough order of magnitude that is a fair 
estimate for galaxies in an intermediate Ldust range (but 
with a large dispersion). Together with the addition of 
a point at z = 6.56 from the upper limit on the gravita- 
tional lens HCM 61 studied by Boone et al. (2007). It is 
important to notice that the dust luminosity of this very 
high redshift object is Ldust < 2.110 11 L Q . So, all the lu- 
minosities lie in the same range - mainly LIRGs. Fig. 13 
seems to confirm the global decreasing trend with the red- 
shift for these galaxies with dust luminosities lower than 
ULIRGs. If we follow the authors who suggest that the 
object is close to a detection, an exponential law is per- 
haps favored. 

8. Testing the MIR Spectral Energy 
Distributions with the /24//15 flux ratio 

The bandpasscs of the AKARI/IRC and Spitzer/MIPS 
filters are larger than the PAH features and do not allow 
us to perform a detailed analysis of the evolution with 
redshift of the MIR SEDs of galaxies. Nevcrthclcs, wc can 
check the consistency of the available libraries with our set 
of data since from z=0 to z=1.2 the IRC 15 /im and MIPS 
24 fim filters encounter the PAH features (cf Fig. 9). Such 
an analysis may also help us to understand the different 
behaviors of the LBG and MIR-selected samples in the 
previous sections where the various calibrations of the 8 
/im flux in bolometric flux are discussed. 

Following Marcillac et al. (2006) and Teplitz et al. 
(2005) we compare the observed 24 to 15 /im flux ratio 
/24//15 with the expected ones calculated to the templates 
presented in the previous section. The predicted fluxes are 
calculated by integrating the templates within the filter 
bandpasses. 

8.1. Variation 0//24//15 with z 

In Fig. 14 we compare the observed variation of /24//15 
with z for both samples (15/im MIR selection and LBGs) 
with the predictions of the differents models presented in 
section 5. The ratio exhibits a large dispersion without 



any significant trend. Error bars are overplotted in the 
sub-panel (c). The templates used for the comparison 
are chosen to sample a large range of luminosity and 
properties. From the CE and LDP libraries, we consider 
four templates corresponding to 10 9 , 10 10 , 10 11 and 10 12 
LQ. Six values of / 6 o//ioo (0,-0.1,-0.2,-0.3,-0.4,-0.5) are 
selected from the DH library in order to cover a large 
range of Ld us t values (from 5 10 8 to 4 10 12 L© according 
to the calibration of Marcillac et al. (2006), /60//100 
increasing monotically with Ldust) 

At z> 0.5 the DH and CE libraries arc able to cover the 
range of observed values although some galaxies at z~l 
exhibit a /24//15 ratio lower than that predicted by the 
models. A substantial fraction of the galaxies detected at 
15 and 24 /zm are LIRGs, this is illustrated in Fig. 14 for 
z>0.9 galaxies (see the next section for the calculation of 
Ldust)- If we restrict the CE templates to Ldust > 10 11 L Q , 
galaxies with /24//15 <~ 1.5 and z>0.9 are not well repre- 
sented by them : templates with a lower Ldust are needed 
to reach such low values of /24//15 which are observed for 
~ 30% of the MIR-selected galaxies with z > 0.9, this issue 
was already underlined by Teplitz et al. (2005) and will 
be re-discussed in the next section. The LDP templates 
predict lower observed ratios which can fit the lowest ob- 
served values but are unable to reproduce the bulk of the 
sample at z>0.5. 

Galaxies at z~l with a large /24//15 ratio can be ex- 
plained cither by a starburst template with a strong con- 
tinuum together with PAHs (CE template with 10 12 L Q , 
DH templates with log(/ 60 //ioo)=0,-0,l also correspond- 
ing to galaxies with Ldust ~ 10 12 L Q ) or by a SED without 
PAHs and a power-law spectrum. This fitting degeneracy 
is also underlined by Marcillac et al. (2006) and Elbaz et 
al. (2005) 

The comparison with the IRS templates implies that the 
spectra of normal star forming galaxies from Smith et al. 
(2007) do not sample the whole set of data and are only 
able to reproduce the lowest observed values whereas the 
starburst template of Brandl et al. (2006) goes through 
the data. The different templates of Spoon et al. (2007) 
exhibit a larger dispersion for the evolution of /24//15 with 
z as expected since they sample extreme classes of objects. 
Templates dominated by the silicate absorption do not re- 
produce at all the general trend of the data: such galaxies 
are rarely found in a MIR selection up to z=l. 

The LBG sample exhibits less dispersed and lower val- 
ues in average of /24//15 than the MIR-selected sam- 
ple. Their behaviour is consistent with that predicted 
for galaxies fainter than 10 12 L© (CE templates) or with 
l°g(/6o//ioo) < — 0.1 (DH templates). The starburst tem- 
plate of Brandl et al. (2006) reproduces well the mean 
behaviour of the LBGs 

To summarize, starburst templates and SED libraries 
for bright IR galaxies (>~ 10 11 L©) are globally consistent 
with most of the /24//15 ratios observed for MIR-selected 
at z>0.5. About ~ 30% of these MIR-selected galaxies 
at z>0.9 exhibit values of /24//15 lower than those for 
starburst and IR bright galaxies and more consistent with 
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Fig. 13. Evolution of L^ ust / Lpjjy ratio for UV-selectcd samples as a function of the redshift from z ~ (Buat et al. 2007) to 
z ~ 1 (Burgarclla et al. 2007) and to z ~ 2 (Rcddy et al. 2006) using the new conversion presented in the text). A point at 
z ~ 3 from Adelberger et al. (2000) and an upper limit from a galaxy at z = 6.56 (Boone et al. 2007). Both an exponentially 
decaying law (dotted line) and a linear regression (dashed line) are a fair representation of the decline of the FUV dust attenuation 
in this redshift range. At z = 0, the red point corresponds to an IR-selected sample and the blue one to an extrapolation from the 
UV-selected sample. 



the ratios found for normal star forming galaxies and tem- 
plates with a moderate luminosity, as well as with the LDP 
library. Nevertheless, most of these galaxies are LIRGs. 
LBGs at z~l exhibit intermediate /24//15 ratios not as 
large as those predicted for bright IR galaxies and galax- 
ies with a strong continuum. They are more typical of 
starburst galaxies as expected for LBGs. 

For z<0.5 galaxies, IRS spectra of normal star-forming 
galaxies (Smith et al. 2007) or SED templates for moder- 
ately bright galaxies (CE, DH, LDP) reproduce the data 
quite well. Strong starbursts or continuum dominated dis- 
tributions (Brandl et al. 2006; Spoon et al. 2007) are 
not consistent with most data. 

8.2. Variation of f 24/ f 15 with L dust 

As underlined above, the various libraries of templates 
can reproduce the full dynamics of the /24//15 ratio but 
not necessarily with the good luminosity. Therefore, we 
have to check their consistency in terms of luminosity. 



Indeed, the CE and LDP templates vary as a function 
of Ijir and the parameter characterizing the DH tem- 
plates, feo/fiooi can be roughly related to an IR lumi- 
nosity (Marcillac et al. 2006). 

As discussed in Sect. 5, the total IR luminosities hdust 
of the observed galaxies are estimated using empirical re- 
lations between hdust and a monochromatic mid-IR lu- 
minosity calibrated on local samples from Takeuchi et 
al. 2005a. Practically, calibrations are available for rest- 
frame data at 12 and 15 ^m. Therefore we will restrict 
our analysis to three redshift bins: 0.2 < z < 0.4 (fluxes 
measured with the IRC 15 fim filter at z = 0.3 correspond 
to 12 //m rest-frame since the isophotal wavelength of the 
L15 filter is 15.6 fim), 0.5 < z < 0.7 and 0.9 < z < 1.1 (fluxes 
measured at z = 0.6 and z = 1 with the MIPS 24 /im filter 
(isophotal wavelength 23.7/^m) correspond respectively to 
15^m and 12/im rest-frame). Using the relations of Chary 
& Elbaz (2001) will give only slightly different luminosi- 
ties which will not have any influence on the following 
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Fig. 14. /24//15 ratio as a function of redshift. Observed ratios are plotted as dots, in panels d,e and f crosses are used for galaxies 
with L^ust > 10 11 Lq and z>0.8, error bars on the measurements are reported in panel c only to avoid crowded plots. The templates 
are overplotted as lines in the different panelsia IRS templates representative of three classes of objects from Spoon ct al. (2007): 
hot dust dominated (1A, solid), PAH dominated (1C (dash), 2C (dot-dash), silicate absorption dominated (3B, dotted); b, IRS 
templates of star-forming galaxies from Smith ct al. (2007): luminosity averaged template (solid) and templates 1 (long dash), 2 
(dotted), 3 (dash) and 4 (dot-dash); c, IRS template for a typical starburst galaxy from Brandl et al. (2006); d, CE templates for 
10 9 (dot-dash), 10 10 (dotted), 10 11 (solid) and 10 12 LQ (dash);e DH templates with / 6 o//ioo=0 (dottcd),-0.1 (dash),-0.2 (solid),-0.3 
(dot-dash) ,-0.4 (lond dash), -0.5 (dash-dot-dot); f: LDP templates, same luminosities and lines as for CE templates (panel d) 
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Fig. 15. /24//15 ratio as a function of L t ; ust for three differents redshift bins. The black dots and error bars are for the MIR-sclectcd 
sample. The LBGs are represented with cyan dots in the third panel, error bars are not represented but are similar to those for the 
IR selected sources. The predictions of the models are overplotted as lines: CE templates (red solid), DH templates (green dash) 
and LDP templates (blue dot-dash) 
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analysis. 

In Fig. 15 are plotted the observed ratios together with 
the predictions of the CE, DH and LDP models for the 
three rcdshift bins. As expected we select brighter galaxies 
when z increases. At low z (0.2 < z < 0.4) the predictions 
of the models are consistent with a large fraction of the 
data but 5 out of the 17 galaxies observed in this red- 
shift bin exhibit ratios larger than those predicted (even 
when error bars are accounted for). Such large values are 
obtained for IRS spectra of strong starburst galaxies as 
can be seen in Fig. 14 (panels a and c). At this redshift 
the observed /24//15 corresponds to /is//i2 m the rest- 
frame of the galaxies: these galaxies are likely to exhibit a 
strong continuum implying a high flux at 18 fim. For the 
medium redshift bin (0.5 < z < 0.7) the models are found 
globally consistent with the data given the error bars on 
the observed ratios. One galaxy is found very discrepant, 
its ratio is only consistent with a template dominate by 
a strong silicate absorption (Fig. 14 (panel a)). At high 
z (0.9 < z < 1.1), both MIR-selected galaxies and LBGs 
are considered since 9 LBGs are detected in our AKARI 
image in the redshift bin 0.9 < z < 1.1. The CE, DH and 
LDP models bracket the observations. Given the large 
dispersions of the observed values the CE and DH tem- 
plates can be considered as globally consistent with the 
data confirming the global reliability of the CE relations 
to estimate Ldust from the monochromatic luminosity at 
8 and 12 /jm rest-frame (see Sect. 6). Half of the sample 
(9/21) correspond to /24//15 < 1-5 whereas the CE and 
DH models predict higher ratios (/24//15 > 1.5) for the 
observed range of luminosity; the LDP models are consis- 
tent with the lowest observed ratios. At z~l the observed 
/24//15 corresponds to /12//8, a s underlined by Elbaz ct 
al. (2005) such low values of this ratio can only be repro- 
duced with a strong PAH emission feature at 7 fim. To 
explain the large variations of this ratio one might invoke 
the relative strengths of the PAH features or the contribu- 
tion of the underlying continuum. On the one hand, the 
galaxies with a low /24//15 (/24//15 <~ 1.5) may have a 
weaker 12 fim PAH emission feature as compared to the 
PAH feature at 7/xm. Indeed Lagache et al. (2004) mod- 
ified the MIR templates in this way and their templates 
can explain the most extreme values of our sample. On 
the other hand, the CE and DH models for bright galaxies 
(as those observed at z~l) exhibit a high continuum which 
contributes to large predicted values of /24//15 without a 
strong variation of the PAH feature at 12 /im; in Fig. 15 
the CE and DH templates of lower luminosity exhibit a 
/24//15 more consistent with the data. In the same way, 
Smith et al. (2007) templates have a rather low contin- 
uum and fit the observed low ratios (Fig. 14 (panel b)). 
Teplitz et al. (2005) also found that a ULIRG spectrum 
is not consistent with the variation of the 24 and 16 /jm 
fluxes observed with Spitzer up to z~1.5 and that local 
star-forming or starburst galaxies are more representative 
of the high z galaxies detected by Spitzer. At z^0.7 Zheng 
et al. (2007) reached a similar conclusion that luminous 
IR galaxies have relatively low temperature, lower than 
that of local galaxies of same luminosity which implies a 



low contribution of the continuum. 

The LBGs at z = 1 also exhibit values of the /24//15 ra- 
tio characteristic of the presence of PAH features in their 
spectra, they appear less dispersed than the IR selected 
galaxies as already underlined. The CE and DH tem- 
plates reproduce rather well the mean values confirming 
the good agreement found in Sect. 5 when Ldust was 
calculated using the CE calibrations at 12 /mum (rest- 
frame). Although the sub-sample considered here is small 
(9 objects) it seems that /24//15 is found higher in LBGs 
with L duat < 10 n L© than for LBGs with L dust > 10 n L©, 
this trend is not seen within the MIR-selected sample and 
most of the MIR-selected galaxies with Ldust < 10 11 have 
a /24//15 ratio lower that those of LBGs. This effect may 
explain the lower L ( j MS t(8/.im) estimated for LBGs than 
for the MIR-selected galaxies for a given L t j, ust (12^m in 
Figs. 11 and 12 (with L dust < 10 n L Q . This different be- 
haviour may be at the origin of the slight differences found 
in the Sect. 6 when the various calibrations of Ldust were 
compared for both samples. 

9. Conclusion 

This paper is based on the observations of the Chandra 
Deep Field South carried out at 15/ito with AKARI/IRC. 
We have applied the standard pipeline to process the 
frame and performed some hotometry to define two galax- 
ies samples with MIR flux densities at A = 15/j.m. A first 
sample is selected at 15/im directly (MIR-selected sam- 
ple) and a second one is based on the mesurement of the 
15/im flux density of an initial sample of 420 LBGs at 
z ~ 1 built by Burgarclla ct al. 2007 with available 24^m 
flux densities from Spitzer/MIPS. 

The main results of this work are : 1) the number 
counts evaluated from the mid infrared-selected are con- 
sistent with previous works from deep imaging at 16^m 
by Teplitz et al. (2005) and Wada et al. (2007). We 
observe a bump around 200/xJ?/ that might be shifted 
to fainter fluxes with respect to other works because of 
Cosmic variance. By dividing the differential counts into 
several sub-samples at different redshifts, we find that the 
bump is mainly due to galaxies at z > 0.4 and at z > 0.8 
for the fainter part of the bump. 2) We have compared 
the dust luminosities estimated from several estimators 
of the monochromatic MIR luminosities at a rest-frame 
wavelength of 8^m with Ldust estimated from MIPS data 
at a rest-frame wavelength of 12/ito assumed to be the 
reference. We find some discrepancies between the esti- 
mates. Bavouzet et al. (2008) and, to some extent, Wu et 
al. (2005) seems to provide dust luminosities consistent 
with the two Ldust estimates from Chary & Elbaz (2001) 
and Takeuchi et al. (2005a). 3) A comparison of the 
/24//15 ratio ratio with several templates and observed 
SEDs is performed that suggest a global consistency of 
the observed ratios albeit there might be an issue with 
the luminosity dependence of this ratio. 4) By making 
use of the best conversion of the MIR luminosity at z = 2, 
we have revisited the evolution of the Ldust/ Luv ratio 
with the redshift and wc find a more drastic decrease of 
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ths raio at z ~ 2 that previously thought by using Reddy 
et al. 2006 conversion. 

10. Acknowledgements 

This work is based on observations with AKARI, 
a JAXA project with the participation of ESA and 
the help of CNES. T.T.T. has been supported by 
Program for Improvement of Research Environment for 
Young Researchers from Special Coordination Funds for 
Promoting Science and Technology, and the Grant-in- 
Aid for the Scientific Research Fund (20740105) commis- 
sioned by the Ministry of Education, Culture, Sports, 
Science and Technology (MEXT) of Japan. A part of this 
work has been done during the stay of D.B and V.B in 
Japan which was supported by Program for Improvement 
of Research Environment for Young Researchers from 
Special Coordination Funds for Promoting Science and 
Technology (MEXT). This paper is the result of a collab- 
oration within the D-SIGALE project. 

References 

Adelberger K. I., Steidel C. C. 2000, ApJ, 544, 218 

Adami, C, Mazure, A., Ilbert, O., Cappi, A., Bottini, D., et 

al. 2005, A&A 443, 805 
Bavouzet, N., Dole, H., Le Floc'h, E., Caputi, K. I., Lagache, 

G. , Kochanek, C. S. 2008, A&A, 479, 83 

Brandl, B.R., Bernard-Salas, J., Spoon, H.W.W. et al. 2006, 
ApJ, 653, 2006 

Boone, F., Schaerer, D., Pello, R., Combes, F., Egami, E. 2007, 

A&A, 475, 513 
Buat, V., Xu, C. K. 1996, A&A, 306, 61 

Buat, V., Iglesias-Paramo, J., Seibert, M., Burgarella, D., 

Chariot, S., Martin, D. C, Xu, C. K., Heckman, T. M., 

Boissier, S., Boselli, A., et al. 2005, ApJ, 619, L51 
Buat, V., Takeuchi, T. T., Iglesias-Paramo, J., Xu, C. K., 

Burgarella, D., Boselli, A., Barlow, T., Bianchi, L., Donas, 

J., Forster, K. 2007, ApJS, 173, 404 
Buat, V., Boissier, S., Burgarella, D., Takeuchi, T. T., Le 

Floc'h, E., Marcillac, D., Huang, J., Nagashima, M., Enoki, 

M. 2008, A&A, 483, 107 
Burgarella, D., Le Floc'h, E., Takeuchi, T. T, Huang, J. S., 

Buat, V., Rieke, G. H., Tyler, K. D. 2007, MNRAS, 380, 

986 

Burgarella, D., Perez-Gonzalez, P. G., Tyler, K. D., Rieke, G. 

H. , Buat, V., Takeuchi, T. T, et al. 2005, A&A, 450, 69 
Caputi, K. I., Lagache, G., Yan, Lin, Dole, H., Bavouzet, N., 

Le Floc'h, E., Choi, P. I., Helou, G., Reddy, N. 2007, ApJ, 
556, 562 

Chary, R., Elbaz, D. 2001, ApJ, 556, 562 

Dale, D.A., Helou, G. 2002, ApJ 576, 159 

Daddi, E., Dickinson, M., Morrison, G., Chary, R., Cimatti, 

A., Elbaz, D., Frayer, D., Renzini, A., Pope, A., Alexander, 

D. M., Bauer, F. E., Giavalisco, M., Huynh, M., Kurk, J., 

Mignoli, M. 2007, ApJ 670, 156 
Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, H., 

Franceschini, A., Fadda, D., Chary, R. R. 2002, A&A, 384, 

848 

Engelbracht, C. W., Rieke, G. H., Gordon, K. D., Smith, J.- 
D. T., Werner, M. W., Moustakas, J., Willmer, C. N. A., 
Vanzi, L. 2008, ApJ 678, 804 



Elbaz, D., Le Floc'h, E., Dole, H., Marcillac, D. 2005, A&A, 
434, LI 

Forster-Schreiber, N., et al. 2004, ApJ 616, 40 
Kawada, M., et al. 2007, PASJ, 59, 389 
Kennicutt R.C. 1998, ARAA 36, 189 

Lagache, G., Dole, H., Puget, J.-L. 2003, MNRAS, 338, 555 
Lagache, G. et al. 2004, ApJS, 154, 112 
Le Floc'h, E. et al. 2005, ApJ 632, 169 

Leitherer, C, Schaerer, D, Goldader, J. D., Delgado, R. M. 

G. , Robert, C, Kune, D. F., de Mello, D. F., Devost, D., 
Heckman, T. M. 1999, ApJS 123, 3 

Lorente, R., Onaka, T., Ita, Y., Ohyama, Y., & 
Pearson, CP. 2007, AKARI IRC Data User Manual, 
ver.1.3, http://www.ir.isas.jaxa.jp/AKARI/Observation/ 

Marcillac, D., Elbaz, D., Chary, R. R., Dickinson, M., 
Galliano, F., Morrison, G. 2006, A&A, 451, 57 

Meurer, G. R., Heckman, T. M., Calzetti, D. 2007, ApJ, 521, 
64 

Murakami, H., et al., 2007, PASJ 59, S369 
Pearson, C. 2005, MNRAS, 358, 1417 

Papovich, C, Rudnick, G., Le Floc'h, E., van Dokkum, P. G., 
Rieke, G. H., Taylor, E. N., Armus, L., Gawiser, E., Huang, 
J., Marcillac, D., Franx, M. 2007, ApJ 668, 45 

Pearson, C, & Rowan-Robinson, M. 1996, MNRAS, 283, 174 

Reddy, N.A., et al. 2006, ApJ 653, 1004 

Rosa-Gonzalez, D., Burgarella, D., Nandra, K., Kunth, D., 
Terlevich, E., Terlevich, R. 2006, MNRAS 379, 357 

Sanders, D. B., & Mirabel, I. F. 1996, ARA &A, 34, 749 

Saunders, W., Sutherland, W. J., Maddox, S. J., Keeble, 
O., Oliver, S. J., Rowan-Robinson, M., McMahon, R. G., 
Efstathiou, G. P., Tadros, H., White, S. D. M., Frenk, C. 
S., Carramiana, A., Hawkins, M. R. S. 2000, AJ 98, 766 

Soifer, B. T., Boehmer, L., Neugebauer, C, Sanders, D. B. 
1989, AJ 98, 766 

Smith, J.D.T., Draine, D.T., Dale, D.A. et al. 2007, ApJ, 656, 
770 

Somerville, R. et al. 2004, ApJ, 600, L171 
Spoon, H.W.W., Marshall, J.A., Houck, J.R. 2007, ApJ, 654, 
L49 

Stoughon, C, et al. 2002, AJ, 123, 485 (erratum 123, 3487) 
Takeuchi, T. T., Ishii, T. T., Hirashita, H., Yoshikawa, K., 

Matsuhara, H., Kawara, K., & Okuda, H. 2001, PASJ, 53, 

37 

Takeuchi, T. T., Buat, V., Iglesias-Paramo, J., Boselli, A., 

Burgarella, D. 2005a, A&A, 432, 423 
Takeuchi, T. T, , Buat, V., Burgarella, D. 2005b, A&A, 440, 

L17 

Teplitz, H. I., Charmandaris, V., Chary, R., Colbert, J. W., 

Armus, L., Weedman, D. 2005, ApJ, 634, 128 
Onaka, T., Matsuhara, H., Wada, T., Fujishiro, N., Fujiwara, 

H. , Ishigaki, M., Ishihara, D., Ita, Y., Kataza, H., Kim, 
W., et al. 2007, PASJ, 59, 401 

Zheng, X. Z., Dole, H., Bell, E. F., Le Floc'h, E. Rieke, G. H., 
Rix, H.-W. Schiminovich, D. 2007, ApJ 670, 301 

Wada, T., Oyabu, S., Ita, Y., Matsuhara, H., Pearson, C. P., 
et al. 2007, PASJ59, 515 

Wolf C, et al. 2004, A&A 421, 913 

Wu H., et al. 2005, ApJ 632, L82 

Xu, C, Lonsdale, C. J., Shupe, D. L., O'Linger, J., & Masci, 
F. 2001, ApJ, 562, 179 



1000 



~i 1 1 — i — i — r 



W 

LO 
i — I 

1—3 

LO 

LO 
* — I 

CO 

lO 

T 1 

Co 



CO 



100 



t 



10 



[] 



~i — i — i — i — i — i 1 1 1 1 — i — i — i — i — i 1 1 



[] 



[] 



[] [] 



1 1 i i i i i 1 1 .i i i i i i i 1 1 1 1 i i i i i i 1 1 1 

1Q -5 1Q -4 1Q -3 1Q -2 

AKARI L15 flux S 15 [Jy] 



